Effective Hamiltonians and dilution effects in kagome and related antiferromagnets 
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What is the zero-temperature ordering pattern of a Heisenberg antiferromagnet with large spin 
length S (and possibly small dilution), on the kagome lattice, or others built from corner-sharing 
triangles and tetrahedra? First, I summarize the uses of effective Hamiltonians to resolve the large 
ground-state degeneracy, leading to long-range order of the usual kind. Secondly, I discuss the 
effects of dilution, in particular to non-frustration of classical ground states, in that every simplex of 
spins is optimally satisfied. Of three explanations for this, the most complete is Moessner-Chalker 
constraint-counting. Quantum zero-point energy may compete with classical exchange energy in a 
diluted system, creating frustration and enabling a spin-glass state. I suggest that the regime of 
over 97% occupation is qualitatively different from the more strongly diluted regime. 



I. INTRODUCTION 

This paper addresses the isotropic antiferromagnet 
with quantum Heisenberg spins Si , on kagome and analo- 
gous lattices ("bisimplex" lattices, to be defined shortly.) 
Everything is restricted to the large- S' and low-T limit, 
where 5* is the spin quantum number and T is the temper- 
ature. Thus, to lowest order we may visualize S^ ~ Ss^, 
where s^ is a classical vector of unit length. 

The theme of Sec. O is the usefulness of effective 
Hamiltonians, in which some degrees of freedom are elim- 
inated, in favor of new terms involving the remaining de- 
grees of freedom. The reader is also reminded that the 
classical picture can be qualitatively wrong, at the tem- 
peratures of interest, which are well below the quantum 
spinwave energies. For example, the quantum pyrochlore 
lattice is seen to be less degenerate than the kagome case, 
in contrast to classical results; and the effective interac- 
tions favoring coUinear/coplanar and other forms of order 
are much more powerful in the quantum case. 

The rest of the paper concerns the effects of dilution. 
Does it produce a spin glass, or generate an effective 
Hamiltonian favoring long-range order, or preserve the 
exceptional degeneracy of the pure lattice? Sec. |l| re- 
views the effect of dilution in ordinary frustrated systems, 
to contrast with its effect in bisimplex lattices (Sec. [rv| ), 
in which the simplex units are all satisfied. Sec. ^ the 
heart of this paper, presents three explanations of the 
simplex satisfaction. But when we admit the full zoo of 
real effects - quantum fluctuations, dilution, unequal ex- 
change constants, anisotropics, external field - it is likely 



that the simplices stop being satisfied (Sec. VC2). 

"Lattices analogous to the kagome" meant, more pre- 
cisely, the bisimplex lattices: those derived from a bipar- 
tite network by placing a spin at each bond midpoint, so 
that each spin belongs to two simplices. (My "simplex" 
is a synonym for "unit" as used by Moessner et al flpl.) 
The coordination number in the network becomes the 
number of corners q of each simplex, which means a sin- 
gle bond, a triangle, or a tetrahedron for g = 2, 3, or 4, 
respectively. Table H lists bisimplex lattices mentioned in 



the literature (there are more) , tagged by the dimension- 
ality and the q, q' values for the simplices on either side 
of a site. 

The kagome and pyrochlore lattices are familiar; the 
connected magnetic sites in SrCr9pGai2-9pOi9 (SCGO) 
form the d — 2 bisimplex lattice that I'll call "kagome 
sandwich" (also known as "pyrochlore slab"). It con- 
sists of two kagome layers connected by a a triangular- 
lattice linking layer ra. The three-dimensional magnetic 
lattice of e.g. gallium gadolinium garnet H] is appropri- 
ately dubbed "hyperkagome" H, since it too consists of 
corner-sharing triangles. The crossed-square lattice is a 
pyrochlore slab normal to {100}, with its top and bottom 
surfaces identified, and serves as a two-dimensional toy 
model for the pyrochlore. (See Fig. 2 of Ref. ||l|(b)) 

The site-percolation threshold Pc of the bisimplex lat- 
tice is obviously the bond-percolation pc of its parent 
bipartite network. It is listed in Table | because one ex- 
pects qualitative changes of behavior at Pc, if the spin 
system orders in any fashion. [|6| The sandwich lattice 
Pc is published here for the first time, to my knowledge. 
(See Appendix 0.) 

The Hamiltonian couples nearest neighbors, 
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with a magnetic field B included. 
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where the exchange constant is J = jS^ > 0, A = 1/2, 
and the total spin of simplex a is 



L„=^i 



(1.3) 



In the kagome-sandwich case, the interlayer (kagome to 
linking layer) coupling will be called J' , but unless explic- 
itly noted, I assume J' = J, the kagome-layer coupling. 



Frustration means generally that, having broken up 
the Hamiltonian into local terms, we cannot simulta- 
neously satisfy all of them. In the present case, each 
term in (1.2) is satisfied (in zero field) whenever Lq, = 0. 



Then each undiluted bisimplex lattice is completely un- 
frustrated from the simplex viewpoint, since it can be 
shown (by example) that every simplex can attain this 
minimum at the same time. Let X denote the classi- 
cal ground state manifold (for B = 0); it is massively 
degenerate on every bisimplex lattice. 



II. EFFECTIVE HAMILTONIANS AND ORDER 
BY DISORDER 



Our goal is to discover how the system breaks symme- 
try and orders (if it does). We assume provisionally that, 
for sufficiently large S, the quantum ground state is one 
of the classical ground states, dressed by some quantum 
zero-point spin fluctuations. (This clearly fails if the fluc- 
tuations about a classical ground state are as big as the 
distance to the next one, or if the system tunnels so freely 
between dressed classical states that the correct wave- 
fmiction is a superposition.) This assumption directly 
gives us the ordering pattern in a simple (e.g. triangular 
lattice) Heisenberg antiferromagnet, where the classical 
ground state is unique (modulo symmetries). 

But commonly, the classical ground states have non- 
trivial degeneracies, so that every ground state has a dif- 
ferent manifold of possible small spin deviations and con- 
sequently a different zero-point energy. Presumably, the 
true quantum ground state should be constructed around 
the particular (ordered!) classical ground state which has 
the lowest zero-point energy. When the thermodynamic 
state of the system has true long-range order in a partic- 
ular pattern, due to this fluctuational energy, we could 
call it "order by disorder" MM . (This usage of the term 
is broader in some ways, and narrower in others, than 
the "local" definition in Ref. Q.) 

To transparently model the degeneracy-breaking, one 
may construct effective Hamiltonians in closed form (via 
intelligently applied perturbation theory). They are in- 
tended, not to provide the accurate energy of a special 
state or two, but as inputs to further modeling, e. g. 
simulations at T > 0, or tunneling calculations M. 



A. Harmonic and higher order effective 
Hamiltonian(s) 

The spin- wave expansion naturally organizes the zero- 
point energy as an expansion in powers of the small 
parameter 1/5', of which the zero term is the 0{jS^) 
classical energy Eq and the first term is the 0{jS) har- 
monic zero-point energy .F(X) = jX^/c^'^fe' including 
all spin-wave frequencies uJk (which depend implicitly on 
X.) This T(X.) is the effective Hamiltonian, defined only 
on states in the ground state manifold. ^ Call the local 
minima of JF() the "favored states" {Y}; these are a dis- 
crete set, (in all cases I know), modulo global rotation 
symmetry. 

Let "ordinary degenerate antiferromagnet" mean one 
for which {X} is a finite-dimensional manifold - clearly 
every ground state is periodic. Examples are the FCC an- 
tiferromagnets pi|jl^ ] or the two-sublattice systems with 
second- neighbor interactions |qJl3,M. (A "highly frus- 
trated" system might be defined as one in which the di- 
mensionality of X is extensive, and two ground states 
may differ only locally.) For an ordinary degenerate sys- 
tem, with isotropic exchange interactions, a crude ap- 
proximation for J^(X) is the biquadratic effective Hamil- 
tonian 
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(2.1) 



This was independently posited phenomenologically [n5J ; 
it can be obtained analytically in a couple of ways using 
a perturbation expansion p2|jlq ] around mean-field the- 
ory. One obtains Kij ~ S^jf'/Shioc, where the local field 
h\oc is 2jS in a bisimplex lattice. Eq. (2.1) correctly tells 
us that the favored states are the collinear ones. In the 



kagome case. (2.1 ) should be replaced by a different func- 



tional form to approximately represent J-"(X), because (i) 
the criterion for favored states is that the spins are copla- 
nar; (ii) the true function JP(Y + (5X) -J^(Y) ~ |X-Y|, 
linear [pp7|] rather than quadratic in the deviations (SX. 
More fundamentally, the harmonic approach fails in 
the kagome case because it does not fully break the de- 
generacy. Every coplanar state has exactly the same 
harmonic-order Hamiltonian, if written using as coordi- 
nates at each site (i) the y (out-of-plane) spin deviation 
and (ii) the rotation angle around the y axis. The source 
of this, mathematically, is that s^ • s^ takes the same value 
(here —1/2) for every nearest-neighbor pair, in every fa- 
vored ground state. Consquently, on a bisimplex lattice 



^Analogously, in the classical model with T/J <C 1 one de- 
fines a harmonic free energy T ^ T <^ Eo ^ J ,hy integrating 
out some degrees of freedom in the partition function. See 
Ref. iioll. 



with only triangles, i.e., the kagome or hyperkagome lat- 
tices, ^(Y) has exactly the same value for every coplanar 
state Y. The number of such states is 0(exp A''), where 
N is the number of unit cells. 

To resolve the surviving degeneracy among favored 
ground states, we need a second effective Hamiltonian 
G(y), obtained from some sort of self-consistent theory 
that takes account of anharmonic spin-wave interactions 
|fi8|-p0[. Since {Y} is discrete, G{Y) is parametrized by 
discrete-spin variables: in the kagome case, the "chiral- 
ities" Tq == ±1 defined from the spins on each plaque- 
tte |l^. Our own approximation Esl^ gave 



(2.2) 



^({'^a}) = -'^JiRa(3)TaT(3, 



a/3 



which has the form of an antiferromagnetic Ising Hamil- 
tonian on the honeycomb lattice, and is defined for every 
favored ground state Y. The energy scale [Q is certainly 
J = 0{jS^^^), i.e. down by only a factor S~^^^ from the 
scale of the harmonic term J^(X). 

In the pyrochlore case, Eq. ( ^.l[ ) is also too crude: 
here the favored ground states Y are the collinear ones, 
but their harmonic energies .?^(Y) are nondegenerate 
since different collinear states have different patterns of 
Si • Sj = ±1. It turns out that the special states that 
minimize J-{) are still infinite in number, but only as 
exp(const L), where L is the system's diameter [ pi[ ; 
the kagome-sandwich lattice seems to behave similarly. 
Thus, in the quantum case, the pyrochlore and sandwich 
lattices show more ordering tendency than the kagome 
lattice (at harmonic order) , whereas in the classical case 
it is the other way around ||l[. 



coupling j, the coefficient in (^^) comes out to Kij = 
jqJS « 7.5K per nearest-neighbor bond. That ought to 
induce a transition to long-range-ordered collinearity at 
a temperature of order lOK), which wouldn't happen at 
any temperature in a classical system [n|. 

III. DILUTION IN ORDINARY CASES 

I now turn to the longer part of this paper: what ef- 
fect(s) does dilution have on a highly frustrated system, 
specifically on a bisimplex-lattice antiferromagnet? Can 
it be represented by an effective Hamiltonian? Or can it 
turn an ordered system into a spin glass? The answers are 
"yes" for the "ordinary" antiferromagnets, which are re- 
viewed in this section (in three regimes) . The answers are 
different for highly frustrated magnets (later sections). 

Starting with this section, I am considering classical 
ground states at T = (unless explicitly noted). 

A. Unfrustrated case 

In an unfrustrated 2-sublattice (Neel) antiferromagnet 
(Fig. H(a)), the ground state at occupied fraction p (on 
the unique extended cluster) has exactly the same spin 
directions as at p = 1, on all the magnetic sites. The total 
magnetization does not cancel exactly, since the moment 
on the even or odd sublattice has 0{\/N) statistical fluc- 
tuations. An observable corollary of this observation is 
that the structure factor at wavevector q. 



'5(^)-i^<iE 



(3.1) 



B. Pitfalls of classical modeling 

Although large S justifies visualizing each spin as a 
fixed-length vector, it does not justify a purely classical 
simulation of the system. The reason is that the interest- 
ing phenomena occur when T <C jS, the scale of JF(X). 
Thus thermal effects are only a small correction to the 
quantum effects. I believe there is an easy fix: a classi- 
cal Monte Carlo simulation using a Hamiltonian which 
includes (an approximation of) .F(X) and ^(Y) ought to 
give valid physical results. 

For a specific example, consider SCGO with S = 3/2 
and a Curie- Weiss constant P4] of 515K, hence j « 80K. 
Then, in a kagome lattice, the effective energy of copla- 
narity would be 0.14jS' per spin along a straight "spin 
fold" 11^] ; if we presume this also applies to the shortest 
rotatable loop of six sites, we get a barrier of lOOK (al- 
ternatively fjS" w 300K from App. B of Ref. §). That 
is vastly larger than the analogous free energy barrier 
in a classical system, at the SCGO freezing temperature 
~ 3.5K. Again, in a pyrochlore lattice with the same 



has the limiting behavior 

lim iS(q) = p{l 

q->0 



p)>0. 



(3.2) 



B. Frustrated case: weak dilution 

Consider next a simple frustrated case, e.g. a Heisen- 
berg triangular antiferromagnet. The pure system has 
a periodic, non-coUinear ground state in which the spins 
differ by 120° angles, which is the best compromise within 
each triangle of spins. The spins in the diluted system no 
longer take the same directions they would in the pure 
system. 

One regime is weak dilution {p close to 1). Consider 
just one nonmagnetic site (Fig. Mb)). Far away from 
this defect, the configuration essentially agrees with a 
pure ground state. But the neighbors of the removed 
spin deviate towards the direction it would have had, and 
their neighbors deviate in turn. Thus the defect creates 
a slowly-varying spin twist with a pseudo-dipolar spatial 



dependence. (It has the angular dependence of a dipole 
field and decays with distance as l/|r|'', where d is the 
spatial dimension.) 

In the ordinary degenerate antiferromagnets, the en- 
ergy reduction due to spin deviations depends on which 
ground state one deviates from. The average of this 
energy, over all the ways to place a low density a;dcf 
of defects, defintes an effective Hamiltonian function 
Hdii(X), proportional to a:def and possessing the full 
symmetry of the ground states X of the pure system. 
In exchange-coupled systems, Hdii(X) favors the least 
coUinear ground states, and thus has an interesting com- 
petition with Hquant §,11 • 

Often, depending on how a defect's local symmetry re- 
lates to that of the ground-state manifold X, each defect 
may prefer a particular ground state, out of the subspace 
{Y} favored by Hdn- In that case, dilution creates an "ef- 
fective random field" MR^ varying in space and coupling 
to the discrete degrees of freedom Y. As in other random- 
field models, disorder wins if it is sufficiently strong or 
the dimension is low enough. The resulting state is a 
spin glass, in the sense that it lacks long-range correla- 
tions and has barriers, but has not been proven to be the 



same phase as a ± J spin glass |23 



C. Frustrated case: strong dilution 

Now return to the same triangular lattice, but with 
strong dilution so that p is just above pc- It is well known 
that, near percolation, the typical connection between 
two sites (if any) is tenuous, and order is propagated 
over one-dimensional chains of sites, which are multiply- 
connected at occasional places. At T = 0, the spin di- 
rections alternate along such a chain, so it constrains the 
relative orientation of the endpoint spins to be parallel or 
antiparallel, depending on whether the number of bonds 
connecting them is even or odd. Order is propagated 
(at T = 0) as if there were a direct bond between the 
endpoints. But if there are two paths in parallel, they 
may disagree on the relation between endpoint spins; the 
smallest example is shown in Fig. yj(c). The ground state 
is a twisted spin-configuration that is found neither in the 
pure lattice nor in a singly-connected chain. In Fig. 0(c), 
the effect is to force an 80° angle between the endpoint 
spin directions. 

The extended connected cluster is an irregular network 
containing loops within loops of this type. It is plausible 
that its global "energy landscape" is like that of a spin 
glass, possessing numerous low-lying, nearly degenerate 
energy minima, separated by energy barriers. 



IV. DILUTION IN A BISIMPLEX LATTICE 



Dilution affects a bisimplex lattice quite differently 
than an ordinary frustrated lattice: simulations find that 
essentially every simplex remains satisfied, i.e. L^ — 0, 
even for strong dilution. (This includes the simplices with 
nonmagnetic sites: they are still simplices, but the num- 
ber of corners q is reduced.) In this section, I will discuss 
the evidence for, and some corollaries of, that fact. 



A. Simulations of 2-layer kagome lattice 

The original evidence was that (in the kagome lattice) 
the local field is exactly 2 on most sites I24|j25ll . I carried 



out more extensive simulations like Ref. |24|, for the di- 
luted kagome-sandwich lattice - relevant to experiments 
on SCGO. (High-quality crystals of that material can be 
grown only with p < 1.) Over 500 independent real- 
izations of the disorder were constructed, for a 10 x 10 
lattice with p = 0.55 (i.e., 385 spins); each realization 
was relaxed from three different random initial configu- 
rations, to a ground state by 250 sweeps, in which each 
site in turn was set to its local-field direction. The pro- 
gram flags all configurations in which |Lq,| > 0.1 on any 
simplex with q > 1. This happened only on "one-eared" 
loops like Fig. Ma). ^ Now, the "one-eared" loop con- 
nects to the rest of the world in (at most) one point. 
Hence it can't induce twisted, frustrated relationships 
among distant spins (fike Fig. ||(c)), and can't be respon- 
sible for spin-glass behavior. Similar results were found 
even when J' ^ J (unequal interlayer and kagome layer 
exchange), as well as in a plain kagome lattice. 



B. "Half-orphan spins" 

Simplex satisfaction has observable consequencies. Ev- 
ery simplex has L^ = - except, of course, that a q = 1 
(one-spin) simplex has JLqI = 1. Let O be the set of 
spins which have g = 1 on one side, with frequency 
Xdci ^ (1 ~ p)^ per unit cell. Also let O' be the spins 
which are completely isolated (q = 1 on both sides) , with 
frequency x'^^^ r^ (1 — p)^ in the sandwich lattice - i.e., 
rare for p > pc- The total magnetization (in units of 
^l = {2hb)S) is 



M. 



9EL- = :^E^»+E 



(4.1) 
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Three exceptions are mentioned at the end of Sec. VB 



The prefactors of 1/2 in (4.1) appear because each spin's 
moment is divided between two simpUces. Schiffer and 
Daruka Q observed a Curie law contribution to the sus- 
ceptibihty, ascribed to "orphan spins" by them (which 
suggests the isolated spins), but more plausibly to the 
spins in O p 27|, which might better be called "half- 
orphan" , since they belong to a simplex on one side but 
are isolated on the other side. If these moments had 



independent directions (but see Sec. IVD), they would 
produce the Curie susceptibility, per cell. 
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(4.2) 



Insofar as the system is built from satisfied simplices, 
its total (classical) magnetization is zero. So, in place of 
(|3.2D, the structure factor (3.2) scales as Xdef ^ (1 —p)- 
Thus even with significant dilution, one expects 5(q -^ 
0) « 0, as seen in SCGO (according to Ref. B). 



The half-orphan site has an free spin of S/2 in (4.1) 
and (^). This is a perfect classical analog of the spin- 
1/2 moment that appears in diluted S = 1 quantum spin 
chains with a valence-bond state pi]. I speculate that 
quantum S — 1 bisimplex lattices indeed have a valence- 
bond state, and S = 1/2 defect moments. 

The simplex-satisfaction concept also encom passes 
bond-randomness effects (see start of Sec. V_C2). The 
ground state of a tetrahedron with one ferromagnetic 
bond still has Lq = so it does not affect the near- 
cancellation of Mtot or the Gauss's law (Sec. IVD). On 



the other hand, a triangle with one ferromagnetic bond 
is unfrustrated, with a configuration like tti- Then 
|Lct| = 1: a ferromagnetic bond produces exactly the 
same sort of paramagnetic defect as a "half-orphan" spin. 



C. Satisfaction with unequal couplings 



What happens to the picture of Eq. (1.2) when J' ^ J 
(coupling to linking layer and/or field B is nonzero? It 
turns out (f^) still works, provided we now take H 






(4.3) 



where Wai = J' / J when a is (before dilution) a tetra- 
hedron and i is the linking-layer spin that caps it, and 
Wai = 1 otherwise. Also, we must replace A in ( |1.2| ) 
by X{a) = J/2J' when a is (before dilution) a triangle. 



and 1 — J/2 J' otherwise (when a is, before dilution, a 
tetrahedron). To satisfy each term, Lq, = A(a)B, which 
immediately implies (in the pure system) 



s. 






B, 



i in kagome layer; 
i in linking layer. 



(4.4) 



Notice that the sum of all linking-layer spins must be ex- 
actly zero either if B — or ii J' — J (even in nonzero 
fiel d). Also, the total magnetization of satisfied simplices 
is (4J), except the coefficients 1/2 must be replaced by 
A(a). Thus, with J' ^ J the net magnetization of the 
satisfied simplices still must be zero. '^ 

There is a problem whenever z is a linking-layer spin 
and a (after dilution) is a (7 = 2 simplex: one must set 
Wai = 1 again to correctly describe the satisfied bond. 
But then, in a magnetic field, the decomposition into 
terms (1.2) breaks down since there is no consistent value 
for A(a). Even for B = 0, if w^i = J' /J still for the other 
simplex that includes site i, the result (4.6) of the next 
section breaks down; in effect, this site enters (4.6) as 
another sort of "point charge" . 



D. Divergence theorem 

I now introduce a sort of "Gauss's Law", which is 
handy for revealing the nonlocal effects of defects. Recall 
that in a "bisimplex lattice", every spin belongs to one 
even and one odd simplex; let (—1)'' = +1 or —1 when /? 
labels an even or an odd simplex, respectively. Mark out 
a domain V containing a subset of the simplices. Define 
a kind of "charge" , 



Q(P)^^(-1)% 



(4.5) 



/3ei? 



The theorem states that, assuming simplex satisfaction. 



^(-i)^«s. = Q{v) = y: (-i)'^'^^'' 



(4.6) 



iedV 



where dV is the set of sites the domain boundary cuts 
through. Also, f3{i) tells which simplex has q = I {oi the 
two containing the half-orphan site i), and S{i) tells which 
simplex is in the interior (of the two containing boundary 
site i). The left-hand side of (4.6) follows since only the 
q — 1 simplices contribute nonzero terms in (4.5); the 
right-hand side follows because every spin in the interior 



^This is contrary to the conclusion of Ref. pi. Their cal- 
culation treated the simplices as independent. That approxi- 
mation can violate important sum rules, e.g. it finds a mean 
spin on the tetrahedron-base which is diflterent from that on 
a triangle, even though these are in fact the same spins. 



of V appears in two terms of ( |4.5| ) with canceling coeffi- 
cients. (Thus, half-orphan spins are the point charges in 
our "Gauss's law", while (— l)''(*^Si plays the role of the 
normal component of the electric field at the surface.) 
This is a generalization of the sum rules of Ref. |l[l(b), 
Sec. Ill B.'* By drawing a succession of nested boundaries 
VI? around a single half-orphan spin s/, one shows that 
neighboring spins Sj have correlations with s/ that alter- 
nate in sign (as speculated by Mendels |Q) and decay 
with distance as l/jr^ — r/|'^~^. 

If we let V include the entire system, there is no bound- 
ary term and the total "charge" (left-hand side of (^)) 
must be zero. This implies that, if there is just one half- 
orphan spin, it is impossible to exactly satisfy every sim- 
plex; if there are just two half-orphan spins, they must be 
exactly parallel or antiparallel (depending on the relative 
parity of their respective q = 1 simplices) . But this law 
has a very large loophole: the violation of L^j = may 
be spread out uniformly over the simplices, such that the 
total energy cost is 0{1/N), which is neghgible in a large 
system. There is a more physical argument why nearby 
"charges" ought, nevertheless, to cancel, as the electro- 
static analogy would suggest. If they don't cancel, spins 
in shells surrounding these "charges" are constrained by 
(4.6) to have a nonzero mean "charge", which would re- 
duce the number of possible states, thus reduce the en- 
tropy, thus increase the free energy at T > 0. I conjecture 
that, in the diluted quantum system, there is an anal- 
ogous effective interaction between nearby half-orphan 
spins, mediated by the harmonic zero-point energy. 



E. NMR experiments 

An NMR experiment measures the distribution of the 
local fields h felt by each NMR nucleus (Ga, in SCGO). 
This is an average of the magnetizations of its neighbor 
(Cr) spins, mi = (s^), where the average is taken over all 
ground states (for a fixed realization), the variance of h 
would scale as Xdof, so the NMR linewidth should scale 
as y/xdei '^ 1 — p, as is seen experimentally 1 27 . 



The experiments can separate the NMR signal from 
the Ga(4f) site, which sees 12 Cr neighbors, including 
three from the linking layer. [p7| Then in the pure system, 
Eq. ( f4.4| ) imphes the mean susceptibility is exactly \i^ jlJ 
per spin, but the mean susceptibilty seen by Ga(4f) is 
Ai'(l/4J-1/8J'). 



V. WHY ARE THE SIMPLICES SATISFIED? 

I now present three different - not exclusive - view- 
points for understanding simplex satisfaction. 

A. Single-impurity explanation 

The original explanation |23] just considered a single 
nonmagnetic impurity in a pure background, as appro- 
priate to the weak dilution regime. Ref. [ p5[ exhibited a 
rearrangement of a few spins (as few as 10) around this 
defect, as in Fig. g, which completely satisfies every sim- 
plex. Farther away, the spin deviations are strictl y zero , 
in contrast to an ordinary frustrated magnet (Sec. Ill B| ) . 
The general method ~ first described for the pyrochlore 
1 rotatable loops, each connecting 



case 



' uses q 

a spin on one of the simplices containing the impurity, 
to a spin on the other affected simplex. ("Rotatable" 
means that, in the pure lattice, the spins on the loop 
can be rigidly rotated together to produce other ground 
states.) The deviations around the impurity remind me 
of the screening around a test charge in a metal, by the 
high density of excitations at zero energy - in the present 
problem, those excitations are the rotatable loops. 

Clearly this picture works for well-separated impuri- 
ties - but that requires 1 — p to be quite small. In Fig. H, 
all 19 spins must point in the pattern shown, modulo 
rotations; this conflicts with the pattern forced by a sec- 
ond impurity that sits anywhere on the two hexagons in 
Fig. y or the eight hexagons surrounding: a single im- 
purity excludes other impurities on 40 other sites. (The 
above counts and Fig. |2| assume a background consist- 
ing of the v3 X v3 state; in other coplanar backgrounds, 
more spins are affected.) 

The single-impurity picture, then, breaks down at 
p « 0.97, where the defect configurations start to over- 
lap. To explain the simulations from this viewpoint, one 
is forced to postulate that the defect spin deviations obey 
a nonlinear superposition principle, as magical as that of 
solitons in certain one-dimensional systems. 



B. Constraint propagation explanation 

This picture is most appropriate to the strong-dilution 
limit near p^j where the connected cluster is tenuous. 



*Note also that, in the pure system (no "charges"), a vector 
potential for the "electric field" can be constructed, which is 
uniquely valued if d = 2; in the kagome case, this is just the 
"spin origami" embedding of Ref. pSJ. 



Along a simple chain of sites, spins alternate, propagat- 
ing order as if there were a direct bond between the end- 
points. Now allow a few spins neighboring this path; 
these decorate the path with triangles which I will call 
"ears" (see Fig. H). Each "ear" removes constraints on 
the spins, since the two spins on the path are now con- 
strained merely to differ by a 120° angle. When a path 
includes two or more "ears" , there is no constraint at all 
between the endpoint spins. 



As in Sec. IIIC, the key question is whether two paths 
which rejoin (forming a loop) might propagate mutually 
exclusive constraints. In contrast to the triangular case 
of Fig. 0(c)), every simple loop has even length. (Derived 
for the kagome case in Ref. Q, Appendix; for a general 
bisimplex lattice, it follows from the bipartiteness of the 
underlying network.) 

All loops with the same number of "ears" are equiv- 
alent, since - by the evenness lemma just stated - they 
have the same number (modulo 2!) of ear-less links, at 
which the spins simply invert. Therefore it suffices to 
study loops of length 6 as in Fig. ||. The one-eared loop 
shown is frustrated: the two spins in the ear should differ 
by 120° due to the triangle, but by 180° due to the five 
simple links connecting them. On the other hand, the 
two-eared loop forces the same 180° angle as a simple 
chain (between the endpoint spins where it connects to 
the rest of the world), and the three-eared loop forces the 
same 120° angles as a giant triangle would. 



The simulations described in Sec. EVA found 38 one- 
ear loops in 500 realizations of 100 cells each a,t p — 
0.55, in full agreement with the predicted frequency 
/[18(l-p)*-H30(l-p)9], per ceh. Three other frustrated 
clusters were found, once each: a one-eared 8-ring, and 
(in two variations) a hexagon pair sharing two ears, with 
one added ear. Each of these objects has (at most) one 
connection to the outside world forming the "ear" , and 
thus has does not propagate frustration on larger scales. 



C. Constraint counting 

The preceding discussions detected no frustration in 
either the strong or weak dilution regime, but are insuffi- 
ciently general. The constraint-counting ( "Maxwellian" ) 
approach of Moessner and Chalker |1|] is, I believe, the 
convincing explanation of simplex satisfaction. ^ The 
basic aim is to compute the dimensionality D of the man- 
ifold of states in which all simplices are satisfied; as long 



as this set is non-null (i.e. £) > 0), it is obviously the 
ground state manifold X. 

Now let Xg be the manifold of "generic" simplex- 
satisfied states (having no linear relationships among the 
spins, apart from the K constraints required to satisfy 
all the simplices). The dimensionality of Xg is 

Dg = F-K. (5.1) 

Here F is the number of degrees of freedom, two per 
spin, so F/N = 2pv per unit cell, where t'sito is the num- 
ber of sites per unit cell, as in Table |. Table H gives 
the number of constraints per simplex, KsiraW), where 
q' is the number of magnetic sites remaining after dilu- 
tion. Naively Ksiin{(l') — 3, for the three components of 
Lq, = (or its generalization to nonzero magnetic field). 
However, in zero field when the simplex is a single bond 
[q' = 2), the constraint is simply "S2 = — Si": attach- 
ing S2 does not change Z?, so the added constraints must 
be two to cancel the added degrees of freedom. Also, 
^sim(O) = -R'sim(l) = 0. The total number of constraints 



K/N - Vt.(g)C>^'(l -p)«-«'x,i,„(g') (5.2) 
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per unit cell, where C', is the combinations of q things 
taken q' at a time, and vi^q) is the number of g-corner 
simplices per cell in the undiluted structure. 

Most relevant to propagating order (or frustration) is 
the generic dimensionality Dgc of the extended connected 
cluster; I approximate this by subtracting from Dg the 
contribution to F by isolated sites. The results are 

Dgc'^'-VN - 6p[l - [l^pf] - 6p2(2 -p), 
Dg,^^"^^ = p[U - 12(1 - pf - 2(1 - pf] - 12p^{p^ - 3p 
DgJ'y^/N - 8p[l - [l-pf] - 12p2(l _p)(2 -p), 

for the three lattices. Countcrintuitively, removing a sin- 
gle site from the kagome lattice leaves Dg unchanged, 
but removing an adjacent pair of sites increases Dg by 1. 
In a general bisimplex lattice, after dilution to the point 
Ptri) at which the average simplex is a triangle, further 
dilution ought to (slightly) increase Dgc- (Note ptii ~ 1, 
0.9, and 0.75 for the kagome, sandwich, and pyrochlore 
lattices respectively.) Indeed, Dgc{p) at first decreases 
rapidly with dilution and and tends to level off below 
Ptii, at roughly Dgc{p) ~ 0.2 on the kagome, ~ 1.6 on 
the sandwich (SCGO) lattice, and ~ 1 on the pyrochlore, 
per unit cell. Since Dgc{p) remains positive, we expect 
simplex-satisfied ground states at any p > pc- 



+ 3), 
(5.3) 



A caution, however, is that "Maxwellian" counting is a 
mean-field theory. A threshold at which Dg (p) vanishes is not 
generally the exact threshold, since portions of the structure 
may be under-constrained while others are over-constrained. 
For related issues in elastic percolation, see Ref. |BI|. 



1. Generic/non-generic state 



crossover; 



Do not be misled by the term "generic" ! The typi- 
cal (physical) ground state is "non-generic" if such states 
have a higher dimensionality than the "generic" ones. 
In particular, in the pure kagome lattice, starting with 
the a/S X V3 coplanar state, one can rotate one of ev- 
ery three six-loops ( "weathervanes" ) by an independent 
angle. This ground-state manifold has dimensionality 
Dn — 1/3 per cell which dominates Dg = 0. Such states 
are non-generic because some spins (e.g. every second 
one in a rotatable loop) have exactly the same directions. 
(Indeed, in ground states relaxed from a random spin 
configuration, second-neighbor spins often point in nearly 
the same direction.) For small dilution, the non-generic 
states still dominate; these are precisely the coplanar 
states with rare impurities discussed in Sec. VA. But 



the spin rearrangement at each impurity (Fig. ph touches, 
and thus immobilizes, 10 six-loops, and the frequency of 
impurities is 3(1 — p) per cell. So, I estimate 

A.(p)/iV«^(l-30[l-p]) (5.4) 

as the dimension per unit cell of the non-generic mani- 
fold. The non-generic manifold loses out to the generic 
one at p « 0.97, in the kagome lattice. By contrast, the 
pure pyrochlore lattice has Dg(l) = 2, while it appears 
-D„(l) — 1 (and presumably Dn{p) also plummets upon 
dilution). Thus, only the generic manifold is relevant in 
the pyrochlore. (essentially this paraphrases the absence 
of (local) "order-by-disorder" [Q.) 



2. Back to frustration 

Up to now, I considered the isotropic classical Heisen- 
berg antiferromagnet, with possible dilution. The degen- 
eracies retained by the classical ground state manifold, 
even under dilution, may be broken by various realistic 
perturbations. These can be handled by the constraint- 
counting framework, as in Table O (with extensions to 
accomodate constraints that involve two simplices). 

Bond disorder has been identified in pyrochlore anti- 
ferromagnets (3|,Q and was modeled theoretically Q. 
It can produce "canted local states" |3y] of deviated spins 
in ordinary Heisenberg antiferromagnets (interacting as 
a spin glass |^3l)- ^^ ^ bisimplex lattice, it has a different 
effect: the constraints -K'sim('?) are greatly increased. As- 
suming a fraction xp of ferromagnetic bonds, inserting 
Ksim{q) from Table || into Eq. ( ^ ) gives 



7-jsand. 
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2(1 
2(1 



ISxi.), 
212;^.) 



(5.5) 



in place of (5.3). Eq. (5.5) predicts that xp ~ 0.05 is 
a threshold, beyond which the system becomes overcon- 
strained and (presumably) spin-glassy. 



Quantum fluctuations also impose collinearity in a 
tetrahedron. For easy-axis exchange anisotropy, rather 
surprisingly the q = 3 ground state manifold's dimen- 
sionality is unchanged |3|]. (The constraint count for 
easy-plane exchange anisotropy is the same as for XY 
spins, see Q). Magnetic field is a special case (already 
covered in Sec. iVB). The spins along a chain when 



B ^ alternate between two possible directions, indicat- 
ing that each extra link contributes only two constraints, 
as accounted in footnote c of Table O. Substitution into 
(^.2[) and (5.1) shows that most of these perturbations 



make the generic states overconstrained, and presumably 
glassy, though (in the undiluted cases) one needs to rule 
out possible non-generic ground states. 



VI. CONCLUSIONS 



Large-S*, non-random antiferromagnet should have pe- 
riodic long-range order at T = 0, even on a highly frus- 
trated lattice. A purely classical picture is invalid, even 
for 5 S> 1, when temperature is far below the spin- wave 
energy jS, but this may be fixed up by using the effective 
Hamiltonians in the simulation (Sec. [IB). 

Dilution does not engender a spin-glass in classical 
ground states, as it does in ordinary frustrated magnets, 
since simplices remain satisfied, as confirmed by simu- 
This was understood most gener- 

The 



lations (Sec. |IVA| ). 

ally from constraint-counting arguments (Sec. 
observed spin-glass state in SCGO might be attributed 
to a competition of dilution with the coplanarity ten- 
dency due to quantum fluctuations. Effective Hamilto- 
nians (Sec. O), serve as a first (perhaps only) way to 
model the zero-point energy contribution in a disordered 
lattice - even if the functional form is not quite right - 
since a more exact treatment would be intractable, in the 



presence of the multiple perturbations of Sec. VC2 



Some physical corollaries are deduced about the struc- 
ture factor, paramagnetic susceptibility, and NMR re- 



sponse, including a "Gauss's Law" rule (Sec. IV D) which 
manifests the long-range effects of the paramagnetic 
("half-orphan") spins. Experiments sensitive to perco- 
lation on SCGO ought to be performed with occupations 
below the threshold Pc ~ 0.5, which is estimated here for 
the first time (Appendix ^) . 

I speculated that - in the plain kagome lattice, at 
least - there is a threshold p, ~ 0.98 separating two 
regimes of dilution: nongeneric coplanar states, peppered 
with defects, at p > p*, but at p < p, generic states, 
in which every trace is lost of the coplanar background 
(Sees. VA and VC). This would imply, of course, that it 



is invalid to extrapolate experimental measurements for 
pG (0.9,0.95) up top = 1. '^ 
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APPENDIX A: PERCOLATION THRESHOLD OF 
KAGOME SANDWICH LATTICE 



A very simple program (no spins are involved) gen- 
erated random configurations with a fixed number of 
diluted sites. Having found the number of sites rrij 
in the j-th connected cluster, and letting mi be the 
largest of them, I computed two percolation quantities: 



the "percolation susceptibility" Xp — ^^^ ^jiH "^j' ^^"^ 
Pp = mi/pN, the fraction in the extensive ("infinite") 
cluster. In the limit N —^ c», one expects Xp ^ \p~Pc\~'^ 
on either side oipc, while Pp ~ |p~Pc|'^ on the high side, 
and Pp = on the low side. 

From the behavior of these two quantities, I estimated 
Pc ~ 0.50(2). To probe the systematic error due to size- 
dependence, several sizes were tried (up to 16 x 16, i.e. 
1792 sites before dilution); however, a genuine scaling fit 
was not carried out. 
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TABLE I. Bisimplex lattices. Number of sites Fsito and of simplices v(q) are per Bravais unit cell. 



Name 


Derived from 


Bravais lattice 


d 


feg') 


i^sitc 


m(3) 


v{A) 


Pc 


kagome 


honeycomb 


triangular 


2 


3,3 


3 


2 





0.65 


garnet 


"gyroid" graph 


bcc 


3 


3,3 


6 


4 





> 0.5 ? 


crossed-square 


square 


square 


2 


4,4 


2 





1 


0.50 


pyrochlore 


diamond 


fee 


3 


4,4 


4 





2 


0.39 


kag. sandwich 


- 


triangular 


2 


3,4 


7 


2 


2 


0.50(2) 



TABLE IL Constraint count Ksi-miq) in g-corner simplices 



Case 


g = 2 


g = 3 


g = 4 


Isotropic Heisenberg spins^ 
The same, with spin waves 
One ferromagnetic bond 
Easy-axis exchange anisotropy 
Magnetic field'' 


2 
2 
2 
4 
3<^ 


3 

3*^ 
4 
4 
3 


3 

5 
6 
5 
3 



"" Also the case J' ^ J 

^ For the coplanarity constraint, add +1 for every site with q ■ 

"^ Also KsiuiiX) ~ 2, since Si || B in that case. 

'^ But subtract 1 for every site with g = 2 on both sides of it. 



3 on both sides, 




(c) 



(a) 



FIG. 1. Diluting ordinary antiferromagnets; diluted sites and their bonds are always 
shown dashed, (a). Unfrustrated antiferromagnet: the total moment here is one down spin, 
the excess of odd sites over even sites (b) . Removal of one spin in the triangular Heisenberg 
antiferromagnet, causing neighboring spins to deviate in the direction of the missing spin. 
Spins on the second ring outwards are given the directions they would have in the pure 
lattice, though in reality a small distortion is found at any radius, (c). A loop of sites with 
an odd number of steps introduces random frustration and spin relationships absent in the 
pure lattice. The upper (lower) path favors the endpoint spins to be parallel (antiparallel). 





FIG. 2. One removed spin in the kagome antiferromagnet, and the 10-site defect loop 
induced around it. Circles with dots or crosses are spins pointing directly out of or into 
the paper; "+" and "— " signs indicate spins with an outward or inward component. The 
surroundings are part of a coplanar \/3 x \/3 ground state. 
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(a). 





FIG. 3. (a,b,c). "One-ear", "two-ear", and "three-ear" loops with their spin configura- 
tions (unique, modulo rotations in spin space) . "A" , "B" , and "C" mark three spin directions 
which differ by 120°. 
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